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Abstract: Magnetic resonance imaging (MRI) is a noninvasive 
method for producing tomographic images of the human body 
and it is most often used for the detection of tumors(lesions) and 
other abnormalities in soft tissues such as the brain. MR Images 
often requires application of noise filtering techniques before 
visual inspection or application of noise sensitive post 
processing methods, so a filtering which is desired to 
significantly decrease image noise and simultaneously to 
preserve some fine details in the image is required .Therefore 
good denoising technique shoud be developed on MRI for 
further processing the Lymph node for automatic segmentation.. 
In this paper Non Linear Anisotropic Diffusion Filtering 
denoising technique is discussed and_ its results were 
compared to the existing wavelet based denoising technique. 


Index Terms: MRI, tumor, brain, segmentation, difussion 
filtering, non linear 


I.INTRODUCTION 


Magnetic resonance imaging technique provoides access 
to impotant anatomical and functional information 
through high speed acquisition.In MRI applications,there 
is an intrinsic trade off between signal to noise ration 
(SNR),contrast to noise ratio (CNR) and resolution. 
Depending on specific diagnostic task,high spatial 
resolution and high contrast may be required,whereas for 
image processing applications,a high SNR is usually 
required ,since most of the post processing segmentaion 
algorithms are very sensitive to noise. Nonlinear 
anisotropic diffusion filters are iterative tunable filters 
introduced by Perona and Malik [1] for the removal of 
background noise in images . They have shown that 
diffusion filters can be used to enhance and detect object 
edges within images . It uses local gradients to control 
the anisotropy of the filter. Since these filters smooth or 
enhance MR images and detect edges they might also be 
used for RF correction and intracranial boundary 
detection in MR images [2]. Perona and Malik formulate 
the anisotropic diffusion filter as a diffusion process that 
encourages intra-region smoothing while inhibiting inter- 
region smoothing. 
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In medical imaging we often face a relatively low SNR 
with good contrast, or a low contrast with good SNR. 
Fortunately human visual system is effective in 
recognizing structures even in the presence of a 
considerable amount of noise .But if the SNR is to small 
or the contrast is too low it becomes very difficult to 
detect anatomical structures [3]. 


Low pass filtering reduces the amplitude of the noise 
fluctuations, but also degrades sharp details such as lines 
or edges .This type of filtering does not respect region 
boundaries or small structures, and the resulting images 
appear blurry and diffused. In case of linear diffusion 
filters the price of eliminating the noise is blurring of 
edges,which leads to the undesirable effects on the 
original image. This undesirable effects can be reduced 
or avoided by the design of non linear filters, the most 
common filtering being median filtering .In median filters 
the edges are retained ,but the filtering results in the loss 
of resolution by suppressing fine details. Recently 
Wavelet based methods are becomong popular in 
preprocessing of Bio-images [5]. 


Spatial filtering technique are applied under following 
assumptions they are (i) The image is supposed to consist 
of many regions in which signal is stationary and ergodic 
in mean and the variance (ii)The image noise is assumed 
to be zero mean and gausian distributed .the main proble 
is to find the stationary region.The problem of finding the 
proper stationary area for local signal estimation is partly 
solved bychoosing filters that are able to distinguish 
homogeneous regions from those with edge regions .This 
anisotropic non linear filter assumes the image noise to 
be gaussian distributed 


Recently developed anisotropic diffusion filtering [1], 
which overcomes the major drawbacks of conventional 
spatial filtering and significantly improves the image 
quality. Anisotropic diffusion filtering is a multistage 
smoothing and edge detection scheme, it can be 
formulated as a diffusion process, and encourages intra 
region smoothing in preference to smoothing across the 
boundaries. In this method the estimation about the local 
image structure is guided by knowledge about the 
statistics of the noise degradation and the edge strengths 


[3 ]. 
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This paper is organized as follows. Section II discuss the 
basic mathematical principles of non-linear anisotropic 
diffusion filters and 1-D anisotropic diffusion filtering. 
Section III briefs the extension of 1-D anisotropic 
diffusion filtering to 2-D anisotropic diffusion filtering. 
Section IV discusses the Proposed method of denoising of 
MR Images’ using Non linear anisotropic diffusion 
filtering and Section V gives the results and discussion. 
Finally Section VI gives the conclusion. 


II. MATHEMATICAL PRINCIPLES OF NONLINEAR 
ANISOTROPIC DIFFUSION FILTER 


Smoothing is formulated as a diffusive process, which is 
suppressed or stopped at boundaries by selecting locally 
adaptive diffusion strengths. In any dimension this 
process can be formulated as the equation given below, 
assuming no sinks or sources[3] 


u(x,t) = div(c(x,t)Vu(x,t)) (1) 


The diffusion strength is controlled by C(x, t) .The vector 


X represents the spatial coordinate(s).The variable t is an 


ordering parameter .The function u(x,t) is taken as 


image intensity I (x, t). 


Two different diffusion function has been proposed, they 
are given below 


ea 2 
an VI(x,t 
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where X is the diffusion constant,it is chosen to preserve 
the edge of the object boundary and to reduce noise 
contribution 


The diffusion function C(x, t) depends on the magnitude 


of the gradient of the image intensity . Its a monotonically 
decreasing function, which mainly diffuses within regions 
and does not affect region boundaries at location of higher 
gradients [3]. 


c(x,t)= f ( VI(x, o) (4) 


If the gradient is large, a discontinuity is assumed and the 
diffusion is halted .The parameter « is chosen according 
to the noise level and the edge strength. To understand the 
relation between the parameter AK and the discontinuity 
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valuev ,it is proposed to define as the product , 
c*WVI called flow .The flow functions related to c, 


and c, are plotted in the Figure 2.1 below for various 


values of K andVI. 


diffusion function 
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Figure 2.1: Flow functions for C, and C, with respect to 


« and VI 


The maximum flow is generated at locations with gradient 
VI equal to «.when decreasing below «_ ,the flow 
reduces to zero, because in homogeneous regions only 
minimal or no flow takes place .Above «x the flow 
function again decreases to zero, halting diffusion at 
locations of higher gradients. A proper choice of the 
diffusion function not only preserves, but also enhances 
the edges .The anisotropic diffusion process can be 
described as a piecewise smoothing of the original signal. 
The assumption of piecewise constant or slowly varying 
intensity is a good estimate to model MR image data 
,which comprise smooth regions separated by 
discontinuities ,representing various tissue categories 
distinguished by different proton densities and relaxation 
properties[4]. 


1-D Anisotropic Diffusion Filtering: 

The intensity change in one iteration step is defined as the 
sum of flow contributions between the neighboring pixel 
intensities. The structure is simulated as a network, 
wherein the centre point of pixel represent nodes and are 
linked together by arcs, whose flow characteristics are 
determined by the conductivity function © [4].The 
network structure for the 1-D filter is given in Figure 2.2 
below. 


I(x- AX) I(x} 
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Figure 2.2: The network structure for 1-D diffusion filters 
[2] 


I(x+ A x) 


The 1-D discrete implementation of the diffusion process 
is defined by the following set of equations as 
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right ~ left 
As the first step the signal flow is calculated between 
neighboring nodes as given above. Then the node 
intensities are updated by the local sum of the flow 
contributions are given as 
0 

I(t+A)=IQ+4 = =[(0+AN*(@ -D) (6) 
The effect of anisotropic filtering in 1-dimension along 
the row is shown in Figure 2.3 and Figure 2.4 which 
shows the smoothing effect (a)before and (b)after 
anisotropic diffusion filtering ,considerable noise 
smoothing effects can be seen without smoothing out the 
edges 
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Figure 2.3: One dimensional signal extracted from the T2 
Image (a) before anisotropic filtering (b) after anisotropic 
filtering to show the effect of edge preserving filtering 
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Figure 2.4: Histogram of (a) T2 image before anisotropic 
filtering (b) after anisotropic filtering 


III. 2D NONLINEAR ANISOTROPIC DIFFUSION 
FILTER 


It’s a simple extension of the 1-D_ discrete 
implementation. The 2-D discrete implementation of the 
diffusion process is is defined by the following set of 
equations as 


£16, t) = div{c(x,t)* grad.I(x,t)]=V" lcs t)*VI(x, ”) 
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The signal flow 
computed as 


between neighboring nodes can be 
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east Dai + ® north = D out (8) 
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The node intensities are updated by the local sum of the 
flow contributions are given as 


I¢+A) =a I=I(Q)+A+*@-O +0 -@D) (9) 


Figure 3.1: The network structure for 2-D diffusion filter 


To obtain better isotropy, flow is also calculated between 
diagonally neighbouring pixels, resulting in an eight way 
connected network. The longer distance between diagonal 


neighbors is taken into account by setting Ad to nD [3]. 


IV. PROPOSED METHOD OF DENOISING OF MR IMAGES 
USING NONLINEAR ANISOTRPIC DIFFUSION FILTERING 


The MRI image to be enhanced is preprocessed before 
applying linear anistropic diffusion filtering for selecting 
proper region of interest. After applying filtering 
technique, it has to be checked for its quality with experet 
opinion. This depicted in Fig 4.1. 


Input MRI Image 


Preprocessi 


Linear Anisotropic 
Diffusion Filtering 


Output 
enhanced 


Fig. 4.1. The Flow chart for the proposed method of 
enhancement 


V. RESULTS AND DISCUSSION 


Retrospetive denoising with a nonlinear technique such 
as anisotropic diffusion filtering has been demonstrated 
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to be an attractive option for improving the SNR of the 
Magnetic Resonance Images. Thus Anisotropic diffusion 
filteringhas proved to be particularly effective in 


prefiltering of MR images before the automatic image 
segmentaion procedure 


Figure 5.1. Images after enhancement (a) Original image 
(b) after histogram equalization (c) after applying 
Modified Retinex algorithm (d) after Contrast limited 
Histogram equalization (e) after applying adaptive 
thresholding (f) after applying unsharp masking 


Figure 5.2 (a) Images after enhancing with symlet wavelet 
by taking ROI (b) Image after 2-D anisotropic diffusion 
filtering 


The visual observation is better in 2-D anisotropic 
diffusion filtering than that of wavelet based enhanced 
image and that the existing algorithms like adaptive 
unsharp masking, contrast limited adaptive histogram 
equalization and other algorithms. This is evident in 
Figure 5.1 and 5.2 


VI. CONCLUSION 


The power of the anisotropic smoothing scheme proposed 
herein lies in the fact that it deals with local estimates of 
underlying image structures, which are highly flexible. 
Discontinuities are preserved and their position is not 
affected. Anisotropic diffusion filtering is an accepted 
filtering technique that is well suited for practical use 
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because of its computational speed and algorithmic 
complexity.Hence the proposed filter can yield highly 
efficient noise reduction and the ability to preserve and 
even enhance important image structures for improved 
visual detection of abnormal tissues in medical image 
diagnostics. 
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